Spectropolarimetry of the 3.4 fim. absorption feature in NGC 1068 
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ABSTRACT 

In order to test the silicate-core/organic-mantle model of galactic interstellar dust, we have 
performed spectropolarimetry of the 3.4 /zm C-H bond stretch that is characteristic of aliphatic 
hydrocarbons, using the nucleus of the Seyfert 2 galaxy, NGC1068, as a bright, dusty background 
source. Polarization calculations show that, if the grains in NGC1068 had the properties assigned 
by the core-mantle model to dust in the galactic diffuse ISM, they would cause a detectable rise 
in polarization over the 3.4 /im feature. No such increase is observed. We discuss modifications 
to the basic core-mantle model, such as changes in grain size or the existence of additional 
non-hydrocarbon aligned grain populations, which could better fit the observational evidence. 
However, we emphasize that the absence of polarization over the 3.4 /im band in NGC1068 — 
and, indeed, in every line of sight examined to date — can be readily explained by a population 
of small, unaligned carbonaceous grains with no physical connection to the silicates. 

Subject headings: galaxies: individual: NGC1068 — galaxies: ISM — galaxies: Seyfert — dust, extinction 
— polarization 

1. Introduction of the origin and evolution of the organic materi- 

als in the ISM, both in our galaxy and in others, 
In our galaxy and no doubt others, stars and ^jj^g Qf gj-eat interest 

planets form from components available in the lo- n ^- i i i- t i ^- i 

, . 1- mi • 1 Particularly puzzlmg are the aliphatic hy- 

cal interstellar medium, i he organic compounds , , , • ,i i-rr Tmv«^ 

.. ur uu drocarbons, whose presence m the dittuse ISM 



observed in interstellar space may therefore be the 
first step towards the complex materials that help 
make planets habitable. A greater understanding 



(DISM) of the Milky Way has been established 
through observations of C-H bond stretches near 
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3.4 /iiTL, seen in absorption throug hout the DISM 
of our ow n galaxy and others ("Sandford et al. 
I991I: [Pendleton et al. 1994; Sandford et al. 199 



Imanishi lioOoHRawhngs et al.ll2003l: iRisaliti et al 



liVl 

al. 



2003l:lDartois et al.ll2004lMason et al.ll2004l) . Two 
very different models exist to explain the pro- 
duction of this material: photoprocessing of the 
icy dust grai n mantles found in de nse molecular 
clouds (e.g. Li fc Greenberd 1997 ). and hydro- 
genation in the DISM of small, amorphous carbon 
particles of the ki nd produced in outfl ows from 



evolved stars (e.g. iMennella et al.l 120021 ). As the 



former process would result in the 3.4 iim band of 
DISM hydrocarbons arising in organic mantles on 
silicate grain cores, it is commonly known as the 
"core-mantle" dust model. 

Both of these models are to some extent con- 
sistent with the observational data. For instance, 
laboratory experiments mimicking both dust for- 
mation pathways have succeeded in producing ma- 
terials with 3.4 /im bands which are a good match 
to th e astronomical feat u re (iGreenberg et al" 



I995I : iMennella erall Il999l l2002h . although the 
critically diagnostic 5-10 /im region reveals a 
much closer match between observations and hy- 
drogenated amorphous carbo n (HAG) materi- 
als (.Chiar et a l. 2000; Pend leton fc AUamandola 
20021 ). The presenc e of the 3.4 iim. ba nd in a pro- 
toplanetary nebula ( Ghiar et al. 1998f l shows that 
hydrocarbon solids can be produced in circumstel- 
lar regions, while the attribution of the 6 /im ex- 
cess absor ption in prot ostars to organic refractory 
material (jGibb fc Wh ittct 2002.) would, if con- 
firmed, indicate that ice processing can also pro- 
duce some fraction of the hydrocarbons. The or- 
ganic OCN~ band at 4.62 /im provid es evidence of 



ice processing in such enviro nments (jDemyk et al 
19981: IPendleton et al.|[l999l ). 



Spectropolarimetry provides opportunities for 
clearly discriminating between these scenarios. In 
the core- mantle model, the silicate core/organic 
mantle grains are responsible for the bulk of the 
visual and IR exti nction and all of the polariza- 
tion in the DISM (|Li fc Greenberelll99lt ). As the 
polarization is caused by dichroic absorption, we 
would therefore expect to observe polarization of 
the continuum, with an increase in polarization 
over absorption fea t ures a s sociated with the grains 



e.g. iHough et all Il988l : IPendleton et all Il990l: 



ization must occur simultaneously over the 3.4 
"mantle" feature, and the 9.7 /im absorption char- 
acteristic of the sil icate core, with rough ly com- 
parable efficiency ( Li fc Greenberel 2002 ). Thi s 
was first investigated by Adamson et al. ( 1999I ). 
who compared the polarization over the 3.4 fxm 
feature towards the Galactic center source, IRS 7, 
with that of the 9.7 fxm band towards the nearby 
Galactic center source, IRS 3. The increase in po- 
larization over the 3.4 /im absorption was found 
to be <0.08%, much less than the 0.4% expected 
on the basis of the silicate feature polarization, in 
apparent conflict with the core-mantle model. 

How ever, as pointed out by 'Li fc Greenberd 
(|2002l ). the silicate feature polarization used in 
that work, while arising toward an object at a 
projected distance of about 0.25pc from that to- 
wards which the 3.4 /im polarization was mea- 



sured ( Geballe et al. 19891 ). does not refer to ex- 



actly the same line of sight, admitting the possi- 
bility that the silicate feature in the diffuse ISM 
towards IRS 7 may be less polarized than ex- 
pected. Studies of additional lines of sight (at 
lower signal-to-noise ratio and/or spectral resolu- 
tion) while reaching si milar conclusions to those of 
Adam son et al. , suffer from similar limita- 



Smith et al.ll2000f ). Specifically, increases in polar- 



tions ( Nagata et al.,.1994: Ishii et a l. 2002). More 
recently, Chiar et al. ( 2006| l measured the polar- 
ization of both the hydrocarbon and silicate fea- 
tures towards the Quintuplet Cluster in the Galac- 
tic Center, finding the hydrocarbons to be signif- 
icantly less polarized than would be expected if 
they existed as mantles on the polarized silicate 
grains. 

To provide another data point against which 
the predictions of the core-mantle grain model can 
be compared, and to extend this study to a new 
and different environment, we have used the bright 
nucleus of the Seyfcrt 2 galaxy NGC1068 as a 
source against which to measure the polarization 
of the 3.4 /im band and the continuum around it. 
Rather than the 10 fim silicate feature, we have 
used the degree of continuum polarization to cal- 
culate the enhancement in polarization that would 
be expected over the 3.4 /im feature from a screen 
of elongated, coated silicate grains which polarize 
by the selective absorption of one plane of the in- 
cident light. In Sj2l we describe the line of sight to- 
wards the nucleus of NGG1068 and its suitability 
for this work. The spectropolarimetric observa- 
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tions and their treatment are outlined in f|3l The 
calculations that we have carried out are discussed 
in [21 and the results presented in ij5l This work 
is then summarized in ^ 

2. Dust and dichroic polarization in NGC 1068 

As the archetypal Seyfert 2 galaxy, NGC 1068 
(d=16 Mpc; Ho = 70km s^^ Mpc'^; 1" = 72 
pc) has been the center of much attention. The 
unified model of active galactic nuclei (AGN) im- 
plies that the galaxy harbors a Seyfert 1 nu- 
cleus, obscured from our point of view by a torus 
of dust and molecular ga s, and 10 /im interfer- 



ometry (jJaffe et al. l2004 as well as the detec- 



tion of broad emission l ines in polarized light 
( Antonucci &: Miller Il985l) lave provided strong 
evidence for the existence of such a torus. Co- 
pious other mid-IR and X-ray data also show 
that the nucleus of NGC 106 8 is obscured by 
large amounts of warm dust (e.g . Bock et al.ll2000t 



Matt et al."2000':'Tom qno et al.ll200ll:lAlloin etal 
2000: .Mason et al. 2006|). While the classificati 



of NGC 1068 as a type 2 object suggests that the 
obscuring torus is oriented roughly edge-on to our 
line of sight, the disk of the galaxy has quite a low 
inclination (i w 29 °: ICarcfa-Gomez et al. 2002 ). 
This means that the line of sight to the center of 
NGC 1068 samples dust local to the active nu- 
cleus, with a negligible contribution from dust in 
the disk of the galaxy. 

The 3.4 /im absorption band has been observed 
in NGC 1068 (iBridger et al.l 11994 IWright et al 



1996: 'i manishi et al.lll997t iMarco k Brooksll200a 
Mason et al. as has silicate absorption at 

9.7 um (iKleinmann et aLlll976l:lRoche et al.lll984l: 



Jaffe et al.ll2004l:lMason et al.ll2006HRhee fc Larkin 



20061 ). The 3.4 fxm band has a very similar pro- 
file to that observed in Galactic lines of sight; 
although the band carrier exists in the central re- 
gion of an active galaxy, it apparently undergoes 
little ex tra processing compared to the Galacti c 
DISM toartois et all l2004t iMason et all |2004 . 



NGC 1068 has a nuclear dust spiral which ap- 
pears physically conne cted to the galaxy disk 
(jPogge fc Martinil l2002h . so it is quite possible 
that the hydrocarbons in the nucleus of NGC 1068 
formed in the disk of the galaxy and then mi- 
grated to the center. We therefore proceed on the 
assumption that the hydrocarbon-containing dust 



in the nucleus of NGC 1068 initially formed by one 
of the mechanisms that has been suggested to be 
responsible for Galactic aliphatic hydrocarbons. 

The nucleus of NGC 1068 is known to be po- 
larized from the UV to the mid-IR, and the mech- 
anisms responsible for this po larization have been 
the su bject of intense study. lAntonucci fc Miller 
(|1985[ ) were among the first to observe the UV 
and visible polarization of NGC 1068, and they 
interpreted the high, wavelength-independent po- 
larization as being due to scattering off very small 



partic les, probably free electrons. iBailey et al 
(jl988l ) later discussed the polarization properties 
of NGC 1068 from the UV out to 10/im, and found 
a twist in position angle from the optical into the 
near-IR. They interpreted this as dust or electron 
scattering being replaced by dichroic absorption, 
the preferential absorption of one component of 
the electric vector, as the most important polariz- 
ing mechanism. Their data also show a change 
of ~ 70° in position angle between the L and 
M bands, consistent with dichroic emission from 
aligned dust grains becoming the dominant polar- 
izing mechanism at longer wavelengths. 

Further evidence for the importance of dichroic 
absorption as a polarizing mechanism in the IR 
comes from imaging polarimetry. After analyzing 
the deviations from the centrosym metric patterns 
expec ted from scattering alone, iLumsden et al.l 
()l999l) concluded that their data require another 
mechanism of constant position angle contribut- 
ing to the polarization. The effect of this mecha- 
nism grows with increasing wavelength, as would 
be expected if absorption were beginning to domi- 
nate over scattering, and was estimated to account 
for > 90% of the K band polarization. Lums- 
den at el. were successful in fitting their observed 
JHK polarized flux points with greybody emission 
from hot {T ^ 1200K) dust, reddened by a A^^ '^^ 
extinction curve and scaled by a Serkowski law 
appropriate for moder ately extinguished Gala ctic 
sources. Furthermore, iPackham et al.l (|l997l) in- 
terpreted the different aperture dependence of the 
J-, H- and K- band polarization as consistent with 
a dichroic contribution in the H- and K-bands. 
With higher-resolutio n HST imaging polarimetry, 
Simpson et ah I (l2002h also found that the K-band 
polarization has a contribution from dichroism, 
but they were able to set tighter limits on the lo- 
cation of the dichroic component, to within 1" of 
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the nucleus. 

This agrees with the more detailed modeling of 
Young et al1(|l995l ). who find that neither electron 
nor dust scattering can account for the near-lR 
polarization of NGC 1068. However, if a dichroic 
component is added, the near-IR polarized flux 
spectrum can be reproduced. In this model there 
is still some contribution from electron scattering 
to the near-IR flux, but this decreases with wave- 
length relative to the dichroic component. By the 
K band, electron scattering contributes perhaps 
10% of the total flux, with dust making up almost 
all of the remainder. At L, the scattering contri- 
bution in their model would be negligible. 

Watanabe et al. ( 20031 ) are also successful in 
modeling the near-IR polarization of NGC 1068 
with polarization from aligned dust grains. In ad- 
dition, they raise the possibihty that the IR po- 
larization could be caused by scattering off large 
grains in the torus, and point out that neither 
the 70 "change in position angle nor the absence 
of a centrosymmetric scattering pattern at longer 
wavelengths necessarily rules out scattering as the 
polarizing mechanism. This has yet to be tested in 
any detail, but the implications that it may have 
for the conclusions of this study are discussed in 

m 

While the near-IR data point to dust absorption 
being the dominant polarizing mechanism in the 
L band, spectropolarimetry of the 9.7 /^m silicate 
feature does not show the pronounced polarization 
excess that might be expected if emission or ab- 
sorption fr om aligned silicates were causing the po- 



larization (jAitken et al.l 119841 ). This observation 



was interpreted as evidence that the mid-IR polar- 
ization arises in emission from non-silicate grains 
or has a nonthermal origin. However, the lack of 
polarization in the silicate feature does not neces- 
sarily rule out the presence of aligned silicate-core 
grains in the nucleus of NGC 1068. In radiative 
transfer calculations of dichroic polarization from 
silicate-containing grain mix tures in dusty disks , 
Efstathiou et al.l (|l997l ) and lAitken et all (|2002l ) 
find numerous configurations in which the polar- 
ization over the feature is in fact quite flat, consis- 



tent with the observations of lAitken et al.l (|l984r ). 



Although it appears possible for radiative trans- 
fer effects to suppress the silicate feature polariza- 
tion, the shorter-wavelength 3.4 /im band should 
be much less affected by such interplay between 



absorption and emission. The ratio of the optical 
depths of the 3.4 and 9.7^m bands in NGC 1068 
supports this suggest ion. In the fi ve Galacti c lines 
of sight examined bv lSandford et al.l ( 19951 ) there 
is a fair degree of correlation between the depths 
of the two features (Tg.7/T3.4 = 13 — 19, with the 
higher values towards the Galactic Centre), but 
in NGC 1068 rg. 77x3.4 « 5, suggesting that the 
3.4 ^m band is indeed less affected by underly- 
ing emission than is the silicate feature. This fur- 
ther suggests that treating the dust producing the 
3.4 iim band as a uniform absorbing screen (an 
assumption implicit in the calculations outlined in 
§11]), while undoubtedly a major simplification of 
the dust geometry and temperature structure in 
this AGN, is still a useful approximation. Overall, 
imaging polarimetry, spectral modeling and the 
near-90 ° rotation in position angle all imply that 
a model of the L-band polarization of NGC 1068 
based on selective dust absorption is a reasonable 
representation of the true situation. 

3. Observations and Data Reduction 

L-band spectropolarimetry of the nucleus of 
NGC1068 was obtained on the nights of 2000 
September 19 and 2000 October 6 using the IR- 
P0L2 spectropolarimetry module and CGS4 on 
the 3.8m UK Infrared Telescope on Mauna Kea, 
Hawaii. The 40 1/mm grating and 0.6"-wide slit 
were used, providing R=1360 at 3.4 /im. The 
weather conditions were good during the second 
night, but some thin cirrus was present on the 
firsfl 

In order to obtain Stokes q and u parameter 
spectra a frame was taken with IRP0L2's half- 
wave plate at each of four positions: 0°, 45°, 22.5°, 
and 67.5°. This cycle was repeated 64 times in all. 
At each of these waveplate angles, ordinary and 
extraordinary beams were extracted from the sky- 
subtracted frames, stacked together, and the final, 
total sp ectra combined u sing the "ratio" method 
(see e.g. lTinbergen 1996") to produce Stokes q and 
u parameter spectra. This method has the advan- 
tage of minimizing the effects of variations in sky 
transmission during the observations. 

Prior to calculation of the polarization, further 



^Spectropolarimetry of the 9.7 /im silicate feature was later 
also attempted, but poor weather meant that none of the 
data obtained were useful. 
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treatment of the raw q and u spectra was nec- 
essary. Firstly, data points between 3.3 and 3.4 
jim (observed) were rejected. This part of the 
spectrum contains deep absorptions from the hy- 
drocarbon cement in IRP0L2's WoUaston prism 
(at 3.35-3.4 /im) and from atmospheric methane 
(3.31-3.33 /Ltm), and these hnes do not ratio out at 
all satisfactorily. 

A second effect which must be dealt with is 
a large-amplitude ripple in the spectrum which 
is thought to be caused by multiple reflections 
in the waveplate. The ripple wa s removed using 



the F FT technique described by lAdamson et al 
Briefly, the ripple causes peaks to ap- 



pear around 0.5 and 1.0 times the Nyquist fre- 
quency in the Fourier transform of the spectrum. 
These peaks were set to zero, then the transform 
was inverted. The instrumental zero point of po- 
larization was determined and removed from the 
NGC1068 q and u spectra using observati ons of 
the u npolarizcd star HD18803 (Clavton fc Martin 
Il98l[ ). and the resulting polarization spectrum 
corrected for the imperfect efficiency of the wave- 
plat^. The position angle of polarization (PA) 
was calibrated using observations of the young 
stellar object, AFGL2519, wh ose PA at 3-4 i^m has 
previously been measured bv lHough et al. (|l989V 
An intensity spectrum was also constructed using 
only a few frames taken near in time to HD18803, 
which was used as a telluric standard. 

The polarization, position angle and total flux 
spectra of NGC1068 are shown in Figure [H t o- 



gether with the polarimetri c resuh s of lBailev et al 
(Il988[ ) and lLebofskv et all (Il978[ ). 



4. Polarization Calculations 

To calculate the L-band polarization produced 
by coated dust grains, we have taken advantage 
of the discret e dipole approximation (PDA) code , 
DDS CAT6.1 (|Drainelll988l : lDraine fc Flatai]lll994 
20041 ). In the DDA the dust grain is represented 
as an array of dipole oscillators on a cubic lattice, 
and absorption, scattering and extinction cross- 
sections are then calculated from the amplitudes 
of the dipoles as they interact with the incident 
electric field. 



The optical constants adopted for the grain 

cores and mantles are those of "astronomical sili^ 

cate" (jPraine fc Leell9 84: Drainc 1985; Wcingart ner fc Draind 
20011 ) and the organ ic refractory material of 



Li &: Greenberd (1997). The refractive index of 



as given by 

www.jach.hawaii.edu/UKIRT/instruments/irpol/CGS4/cgs4pol 



the organic component is derived from the Murchi- 
son meteorite and laboratory residues exposed to 
the solar UV field aboard the EURE CA satel- 
lite (|Greenberg et al.lll995t[Greenberg fc Li 199^). 
While both of those materials show absorption 
bands at 5 -9 /xm that are not seen in th e dif- 
fuse ISM (jPendleton fc Anamandolal l2002h . the 
3.4 /zm bands of both bear a good resemblance to 
the observed feature. 

The existence of polarization caused by dust 
absorption indicates that the dust grains must 
be aspherical, but beyond that their shapes are 
not well-established. We have treated the dust 
grains as oblate spheroids with an axial ratio of 
2:1, the cores and mantles having equal eccen- 
tricities. Based on fits to the Trapez ium 9.7/im 
silicate feature, iDraine fc Led (|l984l) find that 
2:1 oblate spheroids are a reasonable rcprcscnta- 
tion o f dust grains, and Hildebrand fc DragovanI 
(|l995l ) also favor oblate over prolate spheroids. 
They conclude that the axial ratio must be <3:1, 
with a best-fitting value of 1.5:1. There has 
been some dispute over whether oblate or prolate 
spheroids give a better fit to the observational data 
(Greenberg fc Li 1996), but we note that the main 
discrepancy in the polarization predicted by these 
two shapes comes over the 9.7/im silicate feature. 
The differences are small at the wav elengths of in- 
teres t in this work (see Fig. 6 of iDraine fc Lee 
Il984l ). The elongation of the grains will affect the 
efficiency with which they polarize, but has little 
effect oi l the wavelength depe ndence of the polar- 
ization ( Greenberg fc Liiri996l ). 

We have adopted simple picket-fence alignment 
for this work (which is equivalent to the per- 
fect spinning alignment approximation for oblate 
grains), with the grains' short axes orientated per- 
pendicular to the direction of propagation of the 
incident light. This has been shown to be quite ad- 
equate if only the wavelength dependence of the 
polariz ation is required, rather than its absolute 
value ( Ghlewicki fc Greenberg! 19901 ). The func- 
tional form used for the grain size distribution is 
htmf^S MEN a~^-^ power law wi th limits of 0.005 and 
0.25 fim (|Mathis et al.lll977f ). 
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Given this information about the size, shape, 
composition and orientation of the grains, DDSCAT 
calculates a number of quantities. Those relevant 
to this work are 



labs — Cabs/'^O-'^ 

} sea — Csca/'^^ 

^ext — Qabs Qsca 

^pol — Qext: II Qext: 



(1) 

(2) 
(3) 
(4) 



in which Qs are efficiencies and Cs cross-sections 
for absorption, scattering and extinction, and 
Qextj II and Qext,-^ refer to extinction efficien- 
cies fo r the two ort hogonal incident polarization 
states (|Drainelll988l ). 

From Qext and Qpoi, extinction and polariza- 
tion spectra for ensembles of grains can be ob- 
tained in the following manner: 



r(A) = N,r 



p{X) = Ngra 



n{a)C, 



geo^ex 



t{a,\)da (5) 



n(a)C, 



!poi 



geo 



(a, X)da (6) 



where Cgeo is the geometrical cross-section of the 
grain in question, Ngrain the column density of 
grains (with n{a), the number of grains in the 
interval a, a -I- da, appropriately normalized) and 



and a„ 



the limits of the size distribution 



to be considered. 

5. Results and discussion 

We are interested in the magnitude of any in- 
crease in polarization that would be expected over 
the 3.4 /im absorption band, given a certain con- 
tinuum polarization caused by dichroic absorption 
by core-mantle grains. To obtain a crude esti- 
mate of the polarization that might be expected 
over the 3.4 /zm band, we can consider the polar- 
ization observed over the 9.7 fim silicate band in 
galactic lines of sight. For a given grain shape, 
size and composition, changes in abundance and 
degree of alignment will alter the level of con- 
tinuum polarization without affecting the ratio 
of feature to continuum polarization, Pfeat/ Pcont- 
In galactic lines of sight that can be fitted with 
pure absorptive polarization, this ratio is typically 
about 3-4, with a couple of objects reaching val- 
ues of -7 (OMCl BN, SgrA GCS3 IL lSmith et J] 



(|2000l )V Alternatively, the optical constants of 
astrono mical silicate give rise to Pfeat/ Pr.nnt ~ 
12 (Dr aine fc Led Il984l : iHildebrand fc Dragovan 
[l995h . The "excess" polarization over an absorp- 
tion featu re is determined by th e strength of the 
band, and lLi fc Greenberg ( 2002 ) have shown that 
silicates and hydrocarbons polarize to a similar de- 
gree per unit optical depth. Obs erved values of 



tqjJtsa range from about 13-19 (Sandford et al 
19m), implying Pfeat/ Pcont -1.15 - 2.0 in the 3.4 
/zm band. For NGC1068, where the L-band con- 
tinuum polarization is about 2.5%, this translates 
to peak polarizations of 2.9 - 5% across the 3.4 
fim. feature. This suggests that a significant ex- 
cess polarization should be detectable over the 3.4 
/im feature, although the spread in silicate feature 
polarizations an d T3.7/T3.4, and the fact that most 



of the sources in ISmith et al.l (|2000[ ) are dominated 
by molecular cloud material mean that these num- 
bers should be treated with some caution. 

To estimate the likely maximum observed 
polarization change over the 3.4 /im band in 
NGC 1068, the best-fitting linear continuum was 
first found for the polarization data. A model 
polarization spectrum was then created from the 
optical depth spectrum derived from the spec- 
tropolarimetry data, assuming that p cx r. The 
model spectrum was normalized to the continuum 
around 3.55 /im, then scaled by various factors and 
the reduced between it and the binned polar- 
ization spectrum calculated for each scaling. The 
best fit corresponds to a deviation of -0.12±0.25% 
(95% confidence) from the linear continuum; the 
data are in fact consistent with a small decrease in 
polarization over the 3.4 /im feature. The model 
fit giving the upper limit on the change in po- 
larization over the 3.4 /im band (i.e., peaking at 
0.13% above the continuum) is shown in Figure [2] 

The polarization calculated as described in 
§ |4] for a population of corc-mantle grains with 
the MRN size distribution is also shown in Fig. 
[2j As grain size may affect the magnitude of 
the 3.4 /im polarization excess ([J5T]), the po- 
larization from a 0.2 /im grain, close to the 
largest size in the MRN distribution and some- 
what larger than the average < a >= 0.087/im 
of th e finite cylinders in the Li fc Greenberg! 
(|l997l ) model, is also shown. The different 
curves in the figure result from different values 
of fcarb, the fraction of the grain volume that 
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is contained in the mantle. Vari ous figures have 
been suggested for this quantity; Li fc Greenberd 
(|2002l ) co nsider fcarb = 0.2,0.5 0.75 w hile ac- 
cordi ng to Chlewicki fc Greenbere (Il990l). fc, 



0.9. Greenberg fc 

(jl997l ) estimate mantle/core volume ratios of 
about 2:1. As expected, the polarization at 
the peak of the 3.4 fim band increases with the 
amount of carbonaceous material in the grains, 
from approximately 3.0% for the 25% -mantle 
MRN grains, to about 3.6% for the 75%-mantle 
particles. 

For the corc/mantle ratios closest to those most 
commonly quoted in the literature, the 3.4 /im po- 
larization predicted by the model is clearly well in 
excess of the observed polarization. In the case of 
the grains with only 25% of the volume contained 
in the mantle, the polarization at the peak of the 
3.4 /im absorption is 0.29% over the linear contin- 
uum fit mentioned above. This is still somewhat 
above the limit we derive on the observed polar- 
ization; the expected and observed polarizations 
are not reconciled even with only a small fraction 
of each grain being composed of organic refractory 
material. The core-mantle grain model, at least in 
the form proposed for the Galactic diffuse ISM, 
is unlikely to explain the L-band polarization of 
NGC 1068. 

5.1. Variations on the basic core- mantle 
model 

Are there likely to be significant differences in 
the grain population(s) in NGC1068 that might di- 
minish the amount of polarization expected over 
the 3.4 /im band, even while core-mantle grains 
are present? The similarity of the band profile 
in NGC 1068 and the Galactic diffuse ISM sug- 
gests that the composition of the hydrocarbon 



mate r ial is similar in both galaxies (jPartois et al 



2004 iMason et al.l 120041) . but other differences 



may a rise. For instance, as the Li fc Greenberel 
(|l997l ) model proposes that the entire infrared 
polarization arises in the core-mantle grains, we 
have so far also assumed this, but if another, non- 
hydrocarbon-containing aligned grain component 
also contributed to the polarization in NGC 1068 it 
would diminish the size of the rise in polarization 
expected over the 3.4 /zm band while contributing 
to Ay but not T3.4. In NGC1068, such a second 
polarizing grain population could conceivably be. 



for example, bare silicates arising from destruction 
of the mantles on some fraction of the grains. 

Values of extinction to the complex, infrared- 
emitting regions of the nucleus of NGC 1068 may 
not be as straightforward to interpret as the ex- 
tinction through the large column of cold dust to- 
wards the Galactic Center, but estimates range 
from ^15 ( Lumsden et al .' 1999; Watanab e et al. 
^003) to ^40 (lYoungeTal. 1995), disregarding 
values based on T3.4. Given that T3.4 ~ 0.1, this 
implies Av/t^a 150 — 400, compared with galac- 
tic values of ~150 (t owards the galactic center; 
Pendleton et al.l 119941 ) or ^several hundred (to- 



wards field stars at various galactic longitudes; 



Rawlings et al.ll2003[ ). There is therefore no com- 



pelling reason to think that grain components 
other than those proposed for the galactic diffuse 
ISM must exist in NGC 1068, but the wide range 
of estimates of Av/t3,4 for both the galactic and 
extragalactic lines of sight means that this cannot 
be ruled out cither. 

Another issue that could affect the 3.4 /im band 
in both extinction and polarization is that of grain 
size. Increasing grain size tends to increase the 
efficiency of continuum extinction while decreas- 
ing the strength of absorption features, so large 
grains might be able to polarize the L-band con- 
tinuum effectively without producing much excess 
polarization through the 3.4 fim band. For grains 
much larger than about 0.2 /im in radius, scatter- 
ing starts to become important {Qabs/Qsca ~ 5 
for a 0.2 /im core-mantle 2:1 oblate spheroid at 
3.3 fim) and calculations of the polarization from 
such particles in a complex system like NGC1068, 
where the inclination of the dusty torus and mul- 
tiple scattering effects may be critical, are beyond 
the scope of this paper. 

It has been suggested that dust grains in AGN 
may be biased towards larger sizes than in the dif- 
fuse ISM, but the evi dence remains contradictory. 
iLaor fc Draind ( 19931 ) pointed out that large (^10 
/im) dust grains are likely to survive longer than 
small grains close to an AGN and will not produce 
a silicate emission feature. The weakness or ab- 
sence of silicate emission in type 1 AGN prompted 
the inclusion of large grains in so me torus models 
van Bemmel fc DuUemond 20031) . but its recent 



discovery in several quasars may argue against 



large grains ()Hao et al.l 120051 : iSiebenmorgen et al 
I2OO5I ). Flat L-M' colors, peculiar Eb-v/Nh and 
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Av/Nh ratios in Seyfert 2 nuclei and the absence 
of the 2 175 A feature in reddened Seyfert Is have 
all been interpreted as eviden ce for large grains 
(|lmanishill200lt iMaiolino et al.ll2001aibi ). but geo- 
metrical effects may also be able to explain many 

of these observations (jWeingartner fc Murrav 

20021). In the spe cific case of NGC1068 , both 



Young et all ^dM) and IWatanabe et alJ (|2003h 



were able to fit the observed J- to K-band and 
optical polarization with grain sizes no different 
from those thought to exist in the galactic DISM, 
although Watanabe et al. suggest that scatter- 
ing from large grains in the torus might be a vi- 
able alternative. Detailed extinct ion and polariza- 



tion calculations such as those of Watanabe et al 



(|2003l ). if extended to the 3.4 jim band, could 
provide valuable constraints on both the size and 
structure of the grains in AGN tori. 

Finally, we note that if the grain size distribu- 
tion is biased towards larger sizes in NGC1G68, 
this could be through preferential destruction of 
small grains, or large grains may have grown by 
coagulation in this warm, dense environment. If 
the latter, then they may have lost some of their 
former core-mantle nature. It seems likely that 
small hydrocarbon inclusions in large, coagulated 
grains would leave their polarization signature on 
the 3.4 /im band, but again, further calculations 
would be needed to test this. 

5.2. Hydrocarbon formation in the diffuse 
ISM 

The absence of excess polarization over the 3.4 
lim band is naturally explained if the feature arises 
in a population of grains that does not contribute 
significantly to the continuum polarization. Such 
a grain population must be small and/or optically 
isotropic ( small grains be ing harder to align than 
large ones: iLazarianllioOS . and references therein), 
and have no physical connection to the silicate 
grains. Recent laboratory work has provided per- 
suasive evidence that small carbon grains like 
those thought to be ejected from AGB stars can 
be hydrog enated during the late r stages of stellar 
evolution dSchnaiter et al. [11999') and in the dif 
fuse ISM (iMennella et aill9 99: Mu noz Caro et al 



20011 iMennella et all |2002| ). Calculations indi- 



cate that hydrogenation proceeds at a fast enough 
rate in the diffuse ISM to balance the dehydro- 
genation caused by UV photons. Conversely, 



in dense molecular clouds, the reduced abun- 
dance of atomic H and the presence of icy man- 
tles act to prevent rehydrogenation of the grains, 
which can still be efficient ly dehydrogenated by 
photo ns and cosmic rays ( Mennella et al.l l200ll 
20031). T hese and other lines of evidence (e.g. 



Shenov et al.ll2003h imply that most hydrocarbon- 
containing grains are formed in the diffuse ISM 
and that re-formation dominates over the dehy- 
drogenation that subsequently occurs, in agree- 
ment with the non-detection of the 3.4 /im feature 
in dense cloud material. This evolutionary sce- 
nario for the aliphatic hydrocarbons is entirely 
consistent with the lack of a 3.4 /Ltm polarization 
excess in NGC 1068 and seve ral Galactic lines 
of sight (Adamson et al. 19991 : 
Chiar et al.ll2006^ . 



Ishii et all 120021: 



6. Summary 

We have performed L-band spectropolarimetry 
of NGC 1068 and shown that the excess polariza- 
tion over the 3.4 /xm feature is below that which 
would be expected on the basis of the silicate- 
core/organic-mantle grain model as applied to the 
galactic diffuse ISM, consistent with a growing 
body of evidence suggesting that the aliphatic 
hydrocarbons in the general diffuse ISM are not 
formed by processing of the ice mantles that form 
on silicate grains in molecular clouds. The coated 
grain model could still be valid in NGC 1068 
if there also exists an extra, non-hydrocarbon 
aligned grain population, or possibly if the grain 
size distribution is biased to larger sizes than in the 
diffuse ISM of our galaxy (detailed calculations of 
both the continuum and 3.4 /im feature polariza- 
tion from micron-sized dust grains may be a useful 
way of constraining the size and/or composition of 
the carbonaceous grain population in AGN). Al- 
ternatively, reaction of small carbon grains with 
atomic hydrogen in the diffuse ISM would be ex- 
pected to produce a population of small, nonpolar- 
izing hydrocarbon-containing grains which would 
naturally explain the lack of polarization of the 
3.4 fim feature. Such a model is also successful 
in accounting for the non-detection of the 3.4 /im 
band in molecular clouds, which is otherwise diffi- 
cult to explain. 
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Fig. 1. — Top: L-band polarization spectrum of 
NGC 1068. The binned spectrum was obtained by 
calculating the standard error of the q and u spec- 
tra in 0.010 nm. bins and is offset by 3% for clar- 
ity. A lso shown are th e measurements oflBailev et al.l 
(|l988l . triangles) and iLebofskv et all l| 19781. square, 
offset by 0.01 fim. in wavelength). Middle: continuum- 
subtracted flux spectrum, showing the peak of the 
3.4 fim absorption feature. Bottom: posi tion an^le of 
polarizatio n, with the me a surem ents of iBailev et al.l 
(fl98a ) and lLebofskv et all (|l978l ). as above. 
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Fig. 2. — Observed and predicted polarization over 
the 3.4 /xm feature in NGC 1068 for an MRN distri- 
bution of coated silicate particles with different thick- 
nesses of mantle. Also shown are the polarization of 
a 0.2 fim grain {fcarb— 0.25), close to the largest size 
in the MRN distribution, and the limit on the change 
in polarization over the feature based on a model fit 
to the observed polarization (see text). The column 
density of grains (eq. 6) was chosen to reproduce the 
observed polarization around 3.55 /im, a region ex- 
pected to be unaffected by the 3.4 /xm feature. 
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